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ABSTRACT 



Aims. We study the magnetohydrostatic equilibrium of magnetic flux tubes with circular cross sections embedded in a 
magnetic surrounding atmosphere. 

Methods. We solve the static momentum equation in 2.5D to obtain the thermodynamics that are consistent with a 
prescribed velocity and magnetic fields. 

Results. We show that force balance is roughly satisfied if the flux tube's magnetic field is aligned with its axis. 
Equilibrium is guaranteed if this magnetic field possesses a transverse component. Several forms of this transverse field 
are investigated. The resulting magnetic field configurations are critically reviewed in terms of the results from spec- 
tropolarimetric observations. The thermodynamic structure that allows the flux tube to be in mechanical equilibrium is 
also calculated. We show that the inferred pressure, density and temperature stratification reproduce intensity features 
similar to dark core penumbral filaments and penumbral grains. 
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0\ 1. Introduction 

t-H . 

04 ■ The filamentary structure of sunspot penumbrae is often ex- 
C*") ' plained in terms of flux tubes (Bellot Rubio 2003; Solanki 
■s^j- , 2003; Thomas & Weiss 2004). The equilibrium configura- 
f^S ■ tion for such penumbral flux tubes has been studied under 
| the thin-flux tube approximation (Montesinos & Thomas 
<^ • 1997; Schlichenmaier et al. 1998 and references therein). 
^ ' This approximation has the advantage that the govern- 
i-H ! ing equations become a set of one-dimensional equations. 
>^ ■ However, the thin-flux tube approximation has limited ap- 
5h ' plicability in the solar photosphere, since the radius of the 
5^ i penumbral flux tubes is typically comparable to the pres- 
sure scale height. 

In this case, the problem becomes two-dimensional: den- 
sity, pressure and temperature vary over the cross section 
of the flux tube. Spruit & Scharmer (2006) and Scharmer 
& Spruit (2006) have suggested that, in this case, homo- 
geneous flux tubes with a circular cross section cannot be 
in mechanical equilibrium. Such flux tubes are subject to 
stretching along the vertical direction and flattening along 
the horizontal direction that changes their shape and will 
likely destroy them. 

In this paper we investigate this conjecture and take 
a first step towards a more realistic modeling of penum- 
bral flux tubes beyond the thin-flux tube approximation. 
We find analytical solutions for the static momentum equa- 
tion for horizontal flux tubes with a circular cross section, 
embedded in a surrounding atmosphere that harbors a po- 
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tential magnetic field. The basic idea is to set up a generic 
magnetic and velocity field that satisfies certain conditions: 
Maxwell equations, boundary conditions, observations etc. 
These magnetic field and velocity vectors are brought into 
the momentum equation. By assuming force balance, we 
are able to obtain the density and pressure (and therefore 
temperature) distribution inside the flux tube. 



This paper is organized as follows: sections 2 and 3 are 
devoted to describing the basic equations and boundary 
conditions between the flux tube and the external atmo- 
sphere. In Section 4 we assume a particular form for the 
magnetic field in the surrounding atmosphere (potential) 
and, applying the boundary conditions, we study which 
configurations for the flux tube magnetic field are plau- 
sible. A general form for the flux tube magnetic field is 
then obtained. In section 5 we apply this general form for 
the flux tube magnetic field to the static momentum equa- 
tion, developing a treatment to obtain the thermodynamic 
structure. Sections 6 presents particular cases for the flux 
tube magnetic field, and discusses: a) the general appear- 
ance of the magnetic field vector and how it compares with 
results from spectropolarimetric observations; b) the ther- 
modynamic structure inside the flux tube; c) solutions of 
the Radiative Transfer Equation in order to simulate the 
observed continuum intensity emerging from such thermo- 
dynamic stratification. In section 7 we critically review our 
results. Finally, section 8 summarizes the main conclusions 
of our work. 
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2. Basic equations and boundary conditions 

Let us consider the case of a horizontal flux tube that carries 
the Evershed flow and is embedded in a static plasma with 
an inclined magnetic field. The flux tube's cross section is 
circular, with a radius R. For simplicity, its central posi- 
tion in the vertical direction will be taken to be at z = 0. 
The stationary momentum equation can be written (in cgs 
units) as 



p(vV)v = -VP g + -j x B + pg 



(1) 



where p and P g are the density and gas pressure. B, j, v 
and g are the magnetic field, velocity, current and gravity 
vectors respectively. The gravity is taken as g = — ge z with 
g = 2.74 x 10 4 cm s~ 2 . Eq. 1 describes the force balance be- 
tween the inertial force, the pressure gradient, the Lorentz 
force and the gravity force. Let us use cylindrical coordi- 
nates with the axis of symmetry along the tube's axis (x- 
direction; radial direction along the penumbra) . According 
to Borrero et al. (2004, 2005 and 2006a; hereafter referred to 
as papers I, II and III), the properties of the flux tubes and 
the background atmosphere change rather smoothly along 
the radial direction in the penumbra. These changes are 
quantitatively smaller than the variations across the plane 
perpendicular to the tube axis, where large gradients can 
be present at very small scales. We will therefore neglect 
the variation of any quantity along the x coordinate. This 
simplifies our problem to a large extent, since now we con- 
sider only dependencies on the plane containing the cross 
section of the flux tube (plane yz; here we can also con- 
sider polar coordinates). We can write, separately for the 
flux tube interior (index 't', tube) and exterior (index 's', 
surroundings) 



B 



B s = B rs (r, 9)e r + Bg s (r, 9)e e + B xs (r, 6»)e x r > R 
B t = B rt (r, 9)e r + B gt (r, 9)e e + B xt (r, 9)e x r<R 

(2) 



= r> R 
v t = w to e x r < R 



Evershed flow 



(3) 



Note that although there is no dependence on the 
^-coordinate, the magnetic field configuration is 3- 
dimensional. The normal component of the external and 
the internal magnetic field are assumed to vanish at the 
tube's boundary r — R. This is important in order to keep 
a clear distinction between the flux tube and the surround- 
ing atmosphere: 



Bl 



Bt 







(4) 



where the superscript * denotes the location of the tube's 
boundary (R, 9). This boundary condition satisfies the con- 
tinuity of the normal component of the magnetic field vec- 
tor across the interface. Substituting our generic magnetic 
and velocity fields (Eqs. 2 and 3) into Eq. 1 we obtain, 



= -%£ + &oB x 

1 dP, 



e g : 



Jx 



Q0 ~l~ c (jxB r jrBx) 



pg sin a 
— pg cos ( 



These equations are valid for both the flux tube's exterior 
and interior. 9 = lies along y > and grows counter- 
clockwise. We can make use of the relation VxB = ^j, 
in order to obtain 



l dB x 
r 89 



Air 

— Jr 



(6) 



3. Pressure and density balance at the boundary 

Both the radial component of external and internal mag- 
netic field vanish in the vicinity of the flux tube's boundary 
(Eq. 4). This leads to continuity of the normal stress across 
the interface (Kippenhahn & Mollenhoff 1975). Therefore, 
by integrating the Maxwell stress tensor across the flux 
tube's boundary we find total pressure balance (thermal 
plus magnetic) between the flux tube's interior and its sur- 
roundings: 



p* 



P, 



(7) 



Bl 

8tt ' 91 ' 8tt 

Taking derivatives with respect to 9 and regrouping terms 
we obtain: 



dP* 
_ gs 

89 



~W ~ 8^89 [Bet + Bxt ~ Bes B ^ (8) 



Terms in the left hand side of Eq. 8 can be written as a 
function of the density using Eqs. 5 and 6, 



8PI 
89 



j*B*R 



p*Rg cos 9 

1 8Bf 
8^ 89 

Now, by substituting Eq. 9 into 8 we obtain 



p*Rg cos 9 



1 



Pt - Ps = 



8 



S,nRg cos 9 89 



si:} 



(9) 



(10) 



Note that Eqs. 7 and 10 link the gas pressure and density 
of the flux tube and the surrounding atmosphere across 
the interface. They have been obtained with only two as- 
sumptions: (a) the radial component of the flux tube and 
surrounding magnetic field has to vanish at the boundary 
r = R and, (b) velocity and magnetic field vectors are con- 
stant in the direction of the tube's axis. 



4. Potential external field 

Let us consider the case in which the external magnetic field 
is potential: V 2 <I> S = and B s = — V$ s . This problem can 
be solved by variable separation. In addition, the boundary 
conditions given by Eq. 4 must be satisfied 



B* sr = (11) 

Consider (see papers I, II; see also Bellot Rubio et al. 2004; 
Bello Gonzalez et al. 2005) that the external atmosphere is 
nearly at rest, and possesses a magnetic field that is some- 
what inclined with respect to the vertical z-axis. We there- 
fore, assume that B s has the following form far away from 
the tube 



Borrero, J.M.: MHS equilibrium of penumbral flux tubes 



3 



lim B s = B cos 7oe z + B sin 7oe x = B cos 70 sin 6e r 

r— >oo 

+ i?o cos 70 cos 6»e + _B sin7 e x (12) 

where B and 70 are the strength and inclination (with 
respect to the vertical to the solar surface) of the external 
field far away from the flux tube. This yields the following 
solution for B s (r, 9) 



Table 1. Parameters of an exponential fit to the verti- 
cal stratification of density and gas pressure for two atmo- 
spheric models representing the sunspot umbra and quiet 
Sun. P and po refer to the gas pressure and density at 
z = 0. 



Model 


Po 


H p 


po 


H P 




[dyn cm~ ] 


[km] 


[g cm - ' ] 


[km] 


Umbra 


2.76 x 10 b 


84 


1.01 x 10~ b 


92 


Quiet Sun 


1.25 x 10 5 


120 


3.18 x 10~ 7 


130 



R 2 

B s = B sin7 e x + B sin 6>cos7 (l 5- ) e r+ 

R 2 " (13) 

B cos6>cos7o(l + ^2")ee 



4.1. Uniform Internal Field 

For the flux tube wc now assume that both the magnetic 
field and Evcrshcd flow are constant and parallel to its axis: 
B t = B to e x and v t = v to e x . By inserting B t (as well as 
the external B s given by Eq. 13) into Eq. 10 we obtain a 
relation between the external and internal densities: 



Pt~Ps = 



B 2 sin cos 2 70 
nRg 



(14) 



Adopting typical values: B ~ 1500 G, 70 — 60°, R ~ 
100 km, we can estimate the density difference across the 
interface to be 



Pt 



Ps 



7xlO" 7 sin0 [gem -3 ] 



(15) 



The physical significance of this equation can be under- 
stood in terms of the total pressure balance given by Eq. 7 
(see also discussion in Borrero et al. 2006b). At the top 
and bottom of the flux tube, the transverse (azimuthal and 
radial) component of the external field (Eq. 13) vanishes: 
B gs (±TT/2,R) = B rs (±Tr/2,R) = 0. This creates a pres- 
sure imbalance at these points that makes the flux tube 
expand vertically (z-direction) and flatten horizontally (y- 
direction). Eq. 14 shows that, in order to avoid this effect 
and keep the flux tube in mechanical equilibrium, the upper 
portions of the flux tube must be denser that the external 
atmosphere: p t > p s such that it becomes antibuoyant. The 
opposite happens in the lower half of the flux tube: p t < p s . 
According to Eq. 15, the density difference needed is about 
7 x 10~ 7 g cm~ 3 for typical penumbral conditions. 

Another way of looking at the problem is in terms of the 
pressure scale heights. Considering total pressure balance 
and assuming that (both inside and outside the flux tube) 
the gas pressure decays exponentially, it is easy to see that, 
for a homogeneous field inside the flux tube, the fact that 
the transverse component of the external field vanishes at 
the top and bottom of the flux tube, implies a larger pres- 
sure scale height inside the flux tube than outside in the 
lower half of the flux tube, but smaller than the external in 
the upper half. 

Table 1 shows the parameters of an exponential fit to 
the vertical density and gas pressure stratification (p(z) <~ 
Po exp — z/Hp, p(z) ~ po exp — z/H p ) for two atmospheric 
models representing the sunspot umbra (Collados et al. 



1994) and the quiet Sun (Borrero & Bellot Rubio 2002). 
Taking these two models as the limiting cases for the exter- 
nal atmosphere in which the flux tube is embedded, we can 
take po from Table 1 and make po = p* s ■ Note that near the 
umbra, where the density at the continuum level is about 
p* ~ 1 x 10~ 6 g cm~ 3 , the existence of a flux tube with a 
uniform magnetic field is plausible. 

Indeed, very similar solutions have been recently pre- 
sented by Schussler & Vogler (2006) in 3D radiative MHD 
simulations (see their Fig. 3) . In those simulations the intru- 
sion into the external field is caused by a nearly field-free 
gap, whereas in this situation the flux tube is filled with 
a horizontal (pointing radially outwards in the penumbra) 
magnetic field. Note that both situations are equivalent in 
terms of the density balance presented in Eq. 10. 

In the outer penumbra, the right-hand side term in 
Eq. 14 becomes smaller as the external field weakens and 
becomes more inclined: B n — 1000 G, 70 — 65°: 



Pt 



p* + 2x 10" 



sm t 



[gem 



(16) 



Here the external density at the continuum level resem- 
bles more closely that of the quiet Sun: ~ 3 x g 
cm~ 3 . On the one hand, the density perturbation (Eqs. 15- 
16) is smaller than the background density. This makes it 
plausible the existence of a flux tube with a homogeneous 
magnetic field and circular cross section. On the other hand, 
those two values are very similar. This places the current 
model at the limit of its validity. There are however, ad- 
ditional effects that contribute to the equilibrium, such as 
the subadiabaticity of the external atmosphere (Borrero et 
al. 2006b), but this only works after the flux tube has been 
already deformed. Since we are interested in keeping a cir- 
cular cross section, we will investigate other possibilities. 
Those could work alone or in combination with a deformed 
tube embedded in a convectively stable atmosphere. 

4.2. Non-uniform internal Field 

Inspection of Eq. 10 reveals that if the flux tube magnetic 
field possesses an azimuthal component such that 



(17) 



'af(6) +4B 2 cos 2 6» cos 2 70 
the density difference (Eq. 10) that follows is 



Pt 



a 1 dF(0) 
SirRg cos 6 AO 



(18) 
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We therefore have the freedom to choose qQ and J- (6) 
in such a way that the right hand side term of Eq. 18 is 
much smaller than the background density. For instance: 
a < 8-rrRgp* and F(0) = sin 9. 

Considering now VBt = we can write 



8{rB rt ) dB, 



ot 



di 



(19) 



Assuming again separation of variables in the following 
form 



B rt (r,9) =H(0)M{r) 
B 0t {r,6) = g(6)M(r) 



(20) 
(21) 



with the following boundary conditions, that are used in 
order to satisfy Eq. 4 and Eq. 17 



A-t(0) = M(R) = A/"(0) = (22) 
Af(R) = 1 G{9) = ^af{9) + 4B 2 cos 2 9 cos 2 7o (23) 

Using Eqs. (19)-(23) and after some algebra, it can be 
shown that a possible magnetic field configuration is given 

by 



B rt (9,r) = ^—(Tl-r) ■ 1 

V ' 2R 2K ' \jaT{9)+ABl cos 2 0cos 2 7o 

dF(0) 



a — — — — 8Bq sin 9 cos 9 cos 2 70 
dp 



(24) 



B 0t {9, r) (3r - 2R) x /aF(9) + 45 2 cos 2 #cos 2 7o 



B T l —B 



■j:I0 



(25) 
(26) 



Note also that Vt is not necessarily parallel to Bt. This 
implies that the magnetic field will change in time according 
to the induction equation. However, once we consider the 
magnetic and kinematic properties of our model (Eq. 3 and 
24 through 26), it turns out that only the diffusive term of 
the induction equation survives, resulting in magnetic field 
changes over very long time scales (~ days), which does not 
affect the discussion presented here. 

The only restrictions in the choice of a and are: a) 
the flux tube transverse magnetic field (given by Eq. 24 and 
25) must not become singular; b) the total transverse field 
must be much smaller than the field along the tube's axis: 
J Bg t + Bf, t « B xt (in order to be consistent with spec- 
tropolarimetric observations ; see Bellot Rubio et al. 2004) 
; c) p[ (Eq. 18) must not become negative (see Sect. 4.1). 



1 Note that a has dimensions of Gauss 2 , where a/8n is equiv- 
alent to a magnetic pressure: dyn cm -2 



5. Thermodynamic structure 

Once the magnetic field in our 2-D domain is known (Eq. 13 
for the external field and Eqs. (24)-(26) for the flux tube 
field) we can apply the static momentum equation (Eq. 1; 
with a velocity field described by Eq. 3) in order to calculate 
a density and pressure distribution. For convenience we will 
work hereafter in Cartesian coordinates. 

Since the external field is potential, j x B = 0, Eq. 1 
translates into 



dP gs (z) 

dy 



= 



Ps{z) = 



ldP gs (z) 



(J 



_gs\ 

dz 



(27) 
(28) 



Therefore it is sufficient to prescribe the behavior of 
the thermodynamic quantities far away from the flux tube 
y — > 00. In our case we have adopted a temperature T s (z), 
density p s (z) and gas pressure stratification P gs (z) as given 
by the cool umbral model of Collados et al. (1994). In order 
to have a better coverage in deep layers this model has 
been extrapolated downwards up to 1 Mm with typical solar 
stratification resulting from 3D numerical simulations (see 
Cheung 2006; Fig. 3.1 and 3.3). 

For the internal (flux tube) atmosphere we must 
consider that the magnetic field configuration given by 
Eqs. (24)-(26) is no longer potential. In particular there 
is a current oriented along the flux tube axis: 



Jt = Jxte x 



_c9_ 

dr 



(rB, 



01 ) 



dB ri 



(29) 



In this case, considering the horizontal (y-axis) direction 
of static momentum equation (Eq. 5) yields the pressure 
distribution inside the flux tube: 



PgtiViz) = ~ J jxtB zt dy + C(z) 



(30) 



where B zt is obtained from Eqs. 24-25 and j x t from Eq. 29. 
The integration constant C(z) can be evaluated by impos- 
ing total pressure balance at the flux tube boundary. 



B 2 sm 2 -f ~B 2 xt0 -aT(9) 



Stt 



P gt (±VR 2 -z 2 , z) = P gs {z) + 

(31) 

Once the pressure is known for every point inside the 
flux tube, the vertical component of the momentum equa- 
tion yields the density: 



Pt(y,z) 



ixtB 



yt 



BP. 



dz 



(32) 



Given the analytical solution for the magnetic field, 
Eqs. (24)- (26), the gas pressure and density also have ana- 
lytical expressions that can be evaluated through Eqs. (30)- 
(32). These solutions ensure that there is no net force acting 
on the flux tube and it is therefore in mechanical equilib- 
rium. Note that the term involving the velocity field (ad- 
vection term) in Eq. 1 vanishes as long as we assume that 
the velocity is concentrated inside the flux tube and that 
it is parallel to its axis (ir-axis; Eq. 3). Once we have two- 
dimensional (yz) maps of the density and gas pressure, the 
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temperature can be evaluated using the ideal gas equa- 
tion. A variable mean molecular weight is used in order 
to account for the partial ionization of the different atomic 
species. 

6. Results for a = 

6.1. Magnetic field 

The first example we will consider is a = 0. This case is 
particularly interesting because it is valid for any T{9). In 
this case, the magnetic field is as follows: 



400 



Mognetic Field Strength 



R 

B s =B sin7 e x + B sin 6>cos7 (l r) e r+ 



,:2 



(33) 



B cos6>cos7q(1 + ^j)ee 



2r 



(34) 



B t =B xta e x - -j^{R - r)B sin6>cos7oe r 
2r 

+ v?> (3r — 2R)B Q cos 6 cos ~/ a eg 
Rr 

Fig. 1 shows the total magnetic field (top panel) and in- 
clination with respect to the vertical (middle panel) result- 
ing from Eq. 33 and Eq. 34 with the following parameters: 
B = 1700 Gauss, 70 = 65°, B xt0 = 1200 Gauss, R = 125 
km. The inclination with respect to the vertical (z-axis) di- 
rection is defined as: 7 = cos~ 1 (B z /B). Bottom panel in 
Fig. 1 displays the field lines for the transverse field (yz 
plane). Inside the flux tube, the magnetic field is mainly 
align with the x-axis: B xt » \/B~f t + B^ t: but is pointing 
slightly upwards along the sides (7 ~ 80°) and downwards 
in its center: 7 ~ 100°. 

In spite of the simplicity of this model, the overall con- 
figuration displays many similarities with recent investiga- 
tions of the penumbral fine structure. For instance: a: at a 
height corresponding to the axis of the flux tube, the inter- 
nal (flux tube) magnetic field is weaker than the external 
(surrounding field) by an average amount of 600 Gauss (pa- 
pers I, II, III; see also Bellot Rubio et al. 2004); b: the aver- 
age flux tube magnetic field is mostly horizontal 7 ~ 90°, 
while the external field is inclined (next to the flux tube) 
about 50° with respect to the vertical; c: an observer at disk 
center could see an inclination for the magnetic field inside 
the flux tube larger than 90° despite the fact that the flux 
tube axis is perpendicular to the observer. This is due to 
the fact that the transverse field inside the flux tube has a 
negative component B z < in its inner part (see Fig.l bot- 
tom panel). The existence of magnetic flux returning into 
the solar photosphere has appeared in many different works 
(see, e.g. Westendorp Plaza et al. 1997; del Toro Iniesta et 
al. 2001; Schlichcnmaier et al. 2004). In the present model, 
returning flux can be explained by either a flux tube point- 
ing downwards or, even if its axis is aligned with the x-axis 
(as in Fig. 1), as a consequence of the transverse field lines: 
B rt and B et . 
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-200 
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Fig. 1. Color plots of the total magnetic field strength 
(top panel) and magnetic field inclination (middle) for 
the following parameters: B — 1700 Gauss, 70 = 65°, 
B x to = 1200 Gauss, R = 125 km. Note that the strength 
of the magnetic field inside the flux tube along its axis is 
B x to = 1200 Gauss, therefore almost all the magnetic field 
seen inside the flux tube goes along the x axis. Bottom 
panel displays the direction of the transverse field. 



6.2. Thermodynamic structure 

Fig. 2 (top and middle panels) shows surface plots of the 
density and gas pressure excess (with respect to the sur- 
rounding at each height) as obtained from the same mag- 
netic field configuration as in Fig. 1. The vertical scale z 



has been shifted so that the the continuum T5 = 1 level 
(denoted by the middle white dashed line) in the surround- 
ing atmosphere lies at z = 0. The center of the flux tube is 
located at zq — km (same height as the continuum level 
in the surrounding atmosphere). 
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As in Sect. 4.1, the upper part of the flux tube is more 
dense that the surroundings at the same height, while the 
opposite happens in the lower part of the flux tube. This 
effect follows from the requirement of total pressure balance 
at the top and bottom of the flux tube. In addition, accord- 
ing to Eq. 18, there is no density change across the bound- 
ary between the flux tube and the magnetic surrounding: 

Ps= Pf 

The temperature distribution is presented in Fig. 2 (bot- 
tom panel). Note that the combination of low density and 
high gas pressure in the lower portions of the flux tube 
yields large temperatures in this region. In addition, the 
high densities found in the upper half of the flux tube re- 
sult in low temperatures there. Depending upon the choice 
of our model parameters Bq, 70, B xt o, zq and i?, sometimes 
unrealistically high/low temperatures (in the lower/upper 
half) may result. The evolution of any hot plasma element, 
that rises adiabatically with typical entropy values corre- 
sponding to the bottom of the convection zone, limits the 
temperature at the photosphere to roughly T max < 13500 
K. Although we are not aware of a lower limit for the tem- 
perature, we will also consider that any value below 2000 
K is not realistic. 

Fig. 3 presents a parameter study aiming at identifying 
the regions where the resulting temperatures are within the 
allowed range (T max < 13500 and T min > 2000) K. Since 
we have 5 model parameters, we have kept three of them 
fixed while varying the other two: Bq and B xt o (top panel), 
R and 70 (middle panel) and finally R and z (bottom 
panel). In these plots, the solid lines separates the regions 
where T max and T m i n are below/above 13500 K and 2000 K, 
respectively. The shaded areas are the regions where both 
constraints are met. Fig. 3 shows that, for very strong fields 
(B > 1800 and B xW > 1600 G), it is difficult to satisfy the 
constraint of T max < 13500 K. The same happens for very 
vertical external magnetic fields (70 < 60°) and for thin flux 
tubes (R < 50 km). However, this does not mean that the 
model presented here is not valid near the umbra. Indeed, 
Fig. 1 shows a magnetic field configuration that very closely 
resembles that region. 

We have also repeated the calculation of Fig. 3 but using 
as a background atmosphere the mean penumbral model of 
del Toro Iniesta et al. (1994). This model is similar (in its 
density and pressure stratification) to a quiet Sun model. 
Our results indicate that temperatures below 13500 K are 
obtained for Bq and B xt o < 1200 Gauss and R > 100 km. 
This is consistent with our discussion in Sect. 4.1. 



6.3. Emergent intensity: dark core penumbral filaments and 
penumbral grains 

High resolution observations of individual penumbral fil- 
aments in the G-band and continuum intensity at 4877 
A (Scharmer et al. 2002; Rouppe van der Voort et al. 2004; 
Sutterlin et al. 2004) show intensity enhancements at the 
edges of the filament. This is accompanied by a drop in the 
core intensity (Fig. 2 and 3 in Scharmer et al. 2002). The 
contrast varies from image to image but in general they are 
larger in G-band images than at 4877 A . Ever since this 
discovery, this particular feature has been referred to as 
dark core penumbral filaments. Other penumbral features 
reveal penumbral grains composed of what looks like three 
different kernels next to each other (Fig. 6 in Rouppe van 
der Voort et al. 2004). 
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Fig. 2. Top panel: Gas pressure excess P g (y, z) — P gs (oo, z). 
Middle panel: Same for density excess p g (y, z) — p gs (oo 1 z). 
Bottom panel: Temperature distribution T{y, z). Note that 
P gs (oo, z) and p gs (oo, z) decay exponentially with height z. 
This distribution ensures that the flux tube is in mechanical 
equilibrium. Dashed white lines indicate (from lowermost to 
uppermost) the T5 = 10, 1, 0.1 levels. 



Borrero, J.M.: MHS equilibrium of penumbral flux tubes 



7 



7 = 60 [Deg]; R = 100 [Km]; Z = [Km] 




1000 1500 2000 

B [Gouss] 



B = 1500 [Gauss]; B XT0 = 1500 [Gauss]; Z = [Km] 

70 I 1 ' 1 1 




R [Km] 



B = 1500 [Gauss]; B XT0 = 1500 [Gauss]; y = 60 [Deg] 

75 ' ' 7 5 " 1 




R [Km] 



Fig. 3. Shaded areas indicate the range for which the model 
parameters results in minimum temperatures (inside the 
flux tube) larger than 2000 K and maximum temperatures 
smaller than 13500 K. Only two parameters are allowed 
to change at a time: B xt o and Bq (top panels), 70 and R 
(middle panel), z and R (bottom panel). 



We have employed a numerical code to solve the 
Radiative transfer equation (Frutiger 2000) and compute 
the values of the observed continuum intensity (at a ref- 
erence wavelength of 5000 A) for each ray path along the 
z-direction. This returns the continuum intensity as a func- 
tion of the position y across the filament. The result is con- 
volved with an Airy function with a FWHM of 100 km in 
order to account for the limited spatial resolution of the 
observations. 



We have made three different experiments in which con- 
tinuum images at a reference wavelength of 5000 A are 
produced by varying, one at a time, different model pa- 
rameters. Results are presented in Fig. 4: varying zq (top 
panel), R (middle panel), 70 (bottom panel). The values 
for the other parameters are indicated in the correspond- 
ing figures. 

In the first experiment (Fig. 4; top panel) the only pa- 
rameter that changes is zq (central position of the flux tube 
with respect to the r = 1 level of the background atmo- 
sphere). The intensity pattern that results is very similar 
to a dark core penumbral filament. The intensity contrast 
between background and bright feature is larger the higher 
the flux tube is located. Note that the background intensity 
is about 0.12 units of the quiet Sun intensity, resulting from 
the use of the cool umbral model by Collados et al. (1994). 

For the second test (Fig. 4; middle panel) we fix z = 
while the flux tube radius R changes from 75 to 150 km. 
For small flux tube radii we obtain an intensity pattern that 
can be identified with a penumbral grain. With increasing 
diameter the emergent intensity morphology shifts to that 
of a dark core filament. In the third test we explore the ef- 
fects of a varying inclination for the external field: 70 (see 
Fig. 4; bottom panel). In this example, more vertical ex- 
ternal fields yield penumbral grains. When 70 increases, the 
intensity at the center of the flux tube decreases and we 
start to see a dark core filament. These are just examples 
of the kind of intensity features that the model is able to 
produce. In all examples we have been careful to select only 
a set of parameters that keep the minimum/maximum tem- 
peratures inside the flux tube above/below 2000/13500 K, 
respectively. 

Inversions of spectropolarimetric observations (papers 
I,II,III; Bellot Rubio et al. 2004; Riiedi et al. 1998, 1999) in- 
dicate that flux tubes are located high above the continuum 
formation height in the inner penumbra but they are found 
near to it in the outer penumbra0. In addition, the external 
field becomes more horizontal as we move away from the 
umbra. Furthermore (and although not entirely clear based 
on the interpretation of the Net Circular Polarization), it 
appears that the radius of the flux tube increases radially 
in the penumbr£0 (re-direction) . To test what kind of ra- 
dial feature would be obtained under these assumptions, 
we have produced 50 different continuum images where the 
flux tube is initially (inner penumbra) located at zq = 50 
but descends linearly down to zq — in the last image (fur- 
ther from the umbra). At the same time, the inclination 
of the external field changes linearly from 60° in the first 
image to 70° in the last one. The radius, R, was kept con- 
stant at 125 km, however. Results are presented in Fig. 5. 
Similarities with observed features in high-resolution con- 
tinuum images (see, e.g. Fig. 5 in Rouppe van der Voort et 
al. 2004) are remarkable. 

Note that for this test we have assumed that some phys- 
ical parameters depend on the x coordinate. This has no 
particular consequence on our main assumptions as long as 



2 A common interpretation of this result is that the flux tube 
is located at a fairly constant height. However, the height separa- 
tion between the flux tube and the external r = 1 level decreases 
away from the umbra as the Wilson depression is reduced. 

3 Magnetic flux conservation, together with a decreasing mag- 
netic field strength inside the flux tube as a function of x, leads 
to a radially increasing flux tube radius. 
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Fig. 4. Variation of the observed continuum intensity at a 
reference wavelength of 5000 A along the cross section of 
the flux tube. Top panel: all parameters are kept constant 
(see text for details) except for the position of the flux tube 
axis with respect the continuum level in the external atmo- 
sphere: zq. Middle panel: same as before but now the radius 
of the flux tube is the only parameter that changes. Bottom 
panel: same as before but now what varies is the inclination 
of the external field 70. 



the changes along the radial direction in the penumbra are 
small (see Sect. 2). 
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Fig. 5. Simulated intensity map of a penumbral filament. 70 
and z change (along the x-direction) from 60° to 70° and 
from 50 to km. The remaining parameters are: Bq — 1700 
Gauss, Btxo — 1300 Gauss and R = 125 km. 



6.4. Results for a ^ 

For a ^ 0, functions such as F{6) = sm.9 or T(9) = cos 9 
lead to singular B r t(r,9) functions, (see Eq. 24). However 
F{9) = cos 2 9, F(6) = sin 2 9 and T{6) = 1 are plausible. 

We have repeated all the calculations presented in 
Sect. 6 for the case a ^ and T(9) = cos 2 9. In this case, 
a is constrained between ~ [— 10 6 ,10 6 ] to avoid singular- 
ity problems in Eq. 24 and large density perturbations in 
Eq. 18. Results show that the maximum temperature inside 
the flux tube increases for a > 0, but the opposite occurs 
for a < (see Fig. 6). This expands the allowed ranges in 
Fig. 3 for a < as compared to a = 0, and vice-versa for 
a > 0. 



7. Critical comparison with observations 

Although the main intent of this paper is to describe equi- 
librium solutions of a horizontal flux tube embedded in 
a magnetic surrounding atmosphere, it is unavoidable to 
make comparisons between the model predictions and ob- 
servations. In previous sections we have mentioned cases 
where that comparison is positive. In this section we turn 
our attention to aspects where our model faces difficulties. 



7.1. Vertical gradients in the field inclination 

In order to reproduce the Net Circular Polarization (NCP) 
observed in sunspots, strong vertical gradients in the field 
inclination (about 40° in about 100 km: Sanchez Almeida 
& Lites 1992) are required. In papers II and III we used 
a simple (the flux tube cross section was assumed to be 
square) flux tube model where those large gradients in 
the inclination were produced at the lower and upper tube 
boundaries. In this manner, we successfully reproduced the 
observed NCP and the full polarization profiles (see also 
Schlichenmaier et al. 2002; Miiller et al. 2002 and Muller et 
al. 2006). 
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Fig. 6. Vertical stratification of the temperature along the 
center of the flux tube: T(0, z) for different values of a. 
Dotted line: a = —7.5 x 10 5 ; solid line: a — 0; dashed 
line: a = 7.5 x 10 5 . For this example we have adopted: 
B Q = 1700 Gauss, B xt0 = 1200 Gauss, 7o = 65°, R = 125 
km, Zq = km. Note that a < produces a smoother 
temperature stratification, with smaller maximum temper- 
atures. The solid thin line corresponds to the surround- 
ing umbral atmosphere (Collados et al. (1994). Diamonds 
represent the extrapolation towards deep layers starting at 
z ~ —150 km (see Sect. 5). 



A closer look to more realistic (circular cross section) 
flux tube models highlights a lack of strong enough vertical 
gradients: see Fig. 10. Since B z = on top of the flux tube, 
the inclination of the external field there is still 90° . A 40° 
change in the inclination along the vertical direction can 
only be attained if we compare the field inclination inside 
the flux tube with the external one 150 km above the up- 
per boundary (compare two vertical dotted lines in Fig. 10; 
middle panel). However, at this point we start moving out- 
side the region in which the spectral lines are formed. 

A possible solution to the problem is to reduce the ra- 
dius of the flux tube. In this case there would be enough 
vertical length to accommodate the flux tube and at the 
same time to let the external field spread upwards long 
enough to see significant changes in the field inclination. A 
second possibility is to relax the hypothesis of circular cross 
sections. A boundary layer with a cusp shape (Schussler & 
Vogler 2006; Scharmer & Spruit 2006) allows B z > on 
top of the flux tube. This of course makes the external field 
more vertical there and increases the A7 across the upper 
boundary. 

7.2. Lower flux tube boundary 

Perhaps the most serious problem that our model faces 
is the location of the lower boundary of the flux tube. 
According to our calculations and examples shown here (see 
Figs. 2 and 6), the T5 = 1 level is always formed inside the 
flux tube. This is a consequence of the high temperatures 
and large opacities present in the lower half of the flux tube. 
If our model is correct, this would indicate that the lower 
boundary of the flux tube would be invisible to spectropo- 
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Fig. 7. Example of the average properties for the flux tube 
(solid) and surrounding atmosphere (dashed) . From top to 
bottom: magnetic field strength, magnetic field inclination 
and temperature. Note the discrepancies and similarities 
with the model for embedded flux tubes used in Borrero et 
al. (2005; Fig. 4). In the middle panel we indicate the sign 
of dj/dz and dv/dz. The sign of the velocity gradient is 
calculated using Eq. 3: the flux tube carries the Evershed 
flow whereas the surrounding is at rest. 

larimetric observations. However, such observations seem 
to indicate (Fig. 7 in Borrero et al. 2005; see also Fig. 6 in 
Borrero et al. 2006a) that the lower boundary is visible. 

Furthermore, the impossibility of observing the flux 
tube lower boundary would make it difficult to explain 
(without calling upon horizontal gradients) the observed 
change in the sign of d"//dz and dv/dz throughout the 
penumbra (see Fig. 5 in Borrero et al. 2004). Indeed, this 
observational result is what led Bellot Rubio et al. (2002) 
and Mathew et al. (2003) to first postulate the existence 
of flux tubes with two defined boundaries. Let us keep in 
mind that, in the flux tube model (Fig. 10; middle panel), 
we have dj/dz > and dv/dz > in the lower half of 
the flux tube, but d-f/dz < and dv/dz < in the upper 
half. It appears now that realistic flux tube models have 
the region where d^/dz > and dv/dz > beyond the line 
forming region. 

8. Conclusions 

We have studied the equilibrium of horizontal magnetic flux 
tubes (with circular cross sections) embedded in an external 
atmosphere with a more vertical magnetic field (assumed 
to be potential). Our study differs from previous investi- 
gations in the fact that the flux tubes arc not necessarily 
thin. We have found that, for typical penumbral conditions, 
the equilibrium condition is hardly satisifed if the flux tube 
magnetic field is aligned with its axis. This situation is al- 
leviated if the tube's magnetic field possesses a transverse 
component. 

We have obtained the thermodynamical configuration 
resulting from different possibilities for the transverse mag- 
netic field inside the flux tube. We have found that the 
thermodynamical structure of such flux tubes reproduces 
the dark core penumbral filaments and penumbral grains 
observed at high spatial resolution. This happens as a con- 
sequence of the high densities and low temperatures found 
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on the upper part of the flux tube (right above the contin- 
uum t = 1 level formation height). 

The overall magnetic configuration explains many (but 
not all, see Sect. 7), results from spectropolarimctric obser- 
vations of sunspot pcnumbrac. It also explains recent obser- 
vations that find the Evershed effect mainly concentrated 
in the dark lanes of penumbral filaments (Bellot Rubio et 
al. 2005; Rimmele & Marino 2006; Langhans et al. 2007). 
Our Fig. 2 shows that the Evershed effect appears predomi- 
nantly at the flux tube's core. This is due to a larger opacity 
at the tube's center (where the dark core originates) as com- 
pared to its sides. Our model also predicts that, although 
weaker, flows should also appear on the brights sides on the 
filaments. This is in agreement with the findings of Bellot 
Rubio et al. 2005. 

Questions such as the stability and temporal evolution 
of the proposed configurations are beyond the scope of this 
work and need to be addressed by other means. In particu- 
lar our approach neglects the energy equation. It is impor- 
tant also to note that the presence of high density plasma 
region immediately above a low density region makes the 
model prone to Rayleigh- Taylor instabilities. It is unclear 
how the horizontal field inside the flux tube might help to 
stabilize the system (Jun et al. 1995). 
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